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ENERGY  ABSORPTION  MECHANISMS  IN  UNIDIRECTIONAL  COMPOSITES 
SUBJECTED  TO  DYNAMIC  LOADING  EVENTS 


INTRODUCTION 


PURPOSE 

This  report  documents  the  numerical  modeling  that  was  conducted  to  study  the  failure 
modes  in  dynamically  loaded,  unidirectional  (UD)  fibrous  composites.  The  purpose  of  this 
research  was  to  investigate  the  energy  absorption  capabilities  of  admissible  failure  modes  for  a 
cross-ply  UD  composite  laminate  subjected  to  severe  dynamic  loading  events  including  blast, 
shock,  and  ballistic  impact. 

Four  failure  mechanisms  were  considered,  namely  fiber  breakage,  matrix  shearing, 
fiber/matrix  debonding,  and  fiber/matrix  delamination  (that  is,  cohesive  interfacial  failure).  Two 
matrix  materials  were  evaluated,  namely  a  rigid  epoxy  and  a  compliant  elastomer  (Kraton-D). 
The  reinforcing  fiber  material  was  ultrahigh  molecular  weight  polyethylene  (UHMWPE). 


BACKGROUND 

Fiber-reinforced  composites  have  become  an  integral  part  of  commercial,  recreation,  and 
defense  markets.  The  proliferation  of  applications  for  fiber-reinforced  composite  technology  can 
be  in  large  part  attributed  to  tailorability — the  single-most  advantageous  aspect  of  fiber- 
reinforced  composites.  Specifically,  tailorability  provides  the  designer  with  unique  levels  of 
customization,  such  as  selection  of  the  fiber  and  matrix  materials,  lamina  architecture,  and 
laminate  stacking  schedules  and  orientations — all  of  which  facilitate  the  design  of  a  composite 
with  optimal  product  performance  with  minimal  weight. 

The  diverse  requirements  for  various  military/defense  composite  applications 
demonstrate  (1)  how  challenging  it  is  to  design  a  functional  composite  that  offers  optimum 
protection  for  vehicles,  vessels,  buildings,  and  even  personnel  and  (2)  how  critical  tailorability  is 
to  that  design  process.  For  example,  military  land  vehicles  are  often  suited  with  hard  armor 
composites  to  withstand  ballistic  impacts  from  enemy  combatants  and  blasts  from  improvised 
explosive  devices  (IEDs);  marine  vessels,  however,  require  protection  from  blast  overpressures, 
shock,  wave  slap,  electromagnetic  threats,  and  corrosion.  Blast-resistance  composites  and 
fabrics  are  being  deployed  in  building  construction,  temporary  shelters,  and  spall  liners  as 
barriers  for  defeating  fragment  threats.  The  National  Aeronautics  and  Space  Administration 
(NASA)  uses  ballistic  composite  laminates  and  fabrics  for  its  deep-space  habitats  to  protect 
against  micrometeor  orbital  debris  (MMOD),  which  travel  at  hyper  velocities  (>  3  km/sec). 
Personnel  armor,  both  soft  and  hard,  range  from  stab-resistant  vests,  to  vests  designed  to 
withstand  .50  caliber  armor-piercing  ammunition,  to  blast-resistance  diapers  for  vehicle 
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personnel  and  for  demolition  and  mine-clearing  troops.  Cavallaro  ’  ’  provides  an  overview  of 
personnel  annors  and  research  that  has  been  conducted  on  plain- woven  soft  body  armors. 

The  growth  of  composites  in  high-performance  markets  continues  to  outpace  the 
development  of  new  and  improved  physics-based  modeling  methods.  Rather,  current  design 
tools  for  severe,  dynamically  loaded  composite  structures,  which  continue  to  lag  at  the 
constituent  and  interfacial  levels,  lack  the  robustness  required  to  suitably  predict  perfonnance 
prior  to  manufacturing  and  prototype  testing.  Because  this  deficiency  confines  the  development 
cycle  to  the  costly  and  time-excessive  build-test-build  approach,  it  is  imperative  that  research  to 
develop  robust  numerical  finite  element  (FE)  tools  capable  of  capturing  stress  wave  propagations 
and  energy  absorption  mechanisms  in  dynamically  loaded  UD  composites  be  pursued. 

The  tailorability  features  afforded  by  composites  include  hosts  of  different  fiber  materials 
(for  example,  aramids,  polyethylenes,  polyesters,  glasses,  and  ceramics)  and  matrix  materials 
(for  example,  polymers,  ceramics,  and  metals)  from  which  the  designer  can  choose.  Also 
available  are  ranges  of  mechanical  properties  for  the  fiber  and  matrix  materials,  such  as  density, 
strength,  elongation,  modulus,  toughness,  strain  rate  dependencies — all  of  which  have  the 
propensity  to  influence  the  stress  wave  propagations  and  energy  absorption  capacities  at  the 
meso-scale.  Crimp  imbalance,  which  represents  one  architectural  modification  that  can  be 
specified/tailored  in  woven  composites  and  fabrics,  can  be  used  to  selectively  alter  stress  wave 
dynamics.'  Similarly,  lamina  stacking  schedules  (for  example,  0/90,  0/45/90/-45,  0/45/0,  etc.) 
and  thicknesses  4' 5' 6  can  be  changed  to  alter  (1)  how  the  stress  waves  propagate  along  the  in¬ 
plane  and  through-thickness  directions  of  the  laminate  and  (2)  the  extent  that  reflections  occur  at 
the  lamina  boundaries  at  the  macro-scale. 

Batra  and  Hassan  developed  custom  FE  software7  to  subject  LTD  composites  to  dynamic 
loading  events  such  as  blast8  and  shock9  for  various  geometric,  loading,  and  material  properties. 
For  LTD  composites  subjected  to  blast,  it  was  shown  that  approximately  15%  of  the  work  done  by 
external  forces  was  dissipated  through  four  failure  modes  and  the  stacking  sequence  strongly 
influenced  the  energy  dissipation.  The  same  conclusions  were  drawn  for  shock  loading,  with  the 
addition  that  fiber  orientation  influenced  both  the  time  and  location  of  failure  mode  initiation  and 
its  direction  of  propagation.  In  both  models,  delamination  served  as  the  damage  mechanism  that 
absorbed  the  most  energy;  however,  Batra  and  Hassan  pointed  out  that  (1)  the  software 
homogenized  the  damage  energies  so  that  the  energies  associated  with  each  failure  mode  were 
easily  obtained  and  (2)  only  three  internal  damage  variables  were  used  but  four  failure  modes 
were  computed.  The  internal  damage  variables  were  used  to  simulate  matrix  cracking,  fiber 
breakage,  and  fiber/matrix  debonding.  A  failure  envelope  was  used  to  model  the  initiation  of 
delamination  of  the  adjoining  layers.  The  failure  modes  not  considered  were  fiber  pull-out,  fiber 
kinking,  fiber  buckling,  and  matrix  crushing. 

Dolce  et  al.10  developed  a  three-dimensional  (3-D)  FE  model  of  a  carbon  fiber-reinforced 
plastic  composite  plate  subjected  to  blast  loads  from  C4  explosive  charges  along  with 
experimental  testing.  The  model,  which  agreed  reasonably  well  with  experimental  data, 
considered  only  delamination  based  on  the  Dycoss  Discrete  Crack  Model  between  the  layers  in 
the  composite  and  the  Chang-Chang  failure  criteria  for  modeling  in-plane  failure  of  the  UD 
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layers.  All  four  failure  modes,  however,  were  not  taken  into  account  in  the  model,  which  may  be 
one  explanation  for  the  difference  in  matching  numerical  results  with  experimental  results. 

Will  et  al.11  investigated  the  effects  of  laminate  stacking  of  carbon  fiber-reinforced 
polymers  subjected  to  projectile  impact.  Their  research  revealed  that,  in  high-speed  impacts,  a 
structure  deformed  locally  and  little  energy  was  used  to  defonn  fibers  and  the  structure. 
Furthermore,  a  significant  amount  of  energy  was  dissipated  in  delamination,  debonding,  and 
fiber  pull-out. 

7  12  13  14 

Other  researchers  ’  ’  ’  developed  cohesive  elements  to  model  damage  initiation  and 
damage  evolution  in  the  fiber/matrix  interfacial  zone.  Karahan15  compared  the  ballistic 
performance  and  energy  absorption  in  woven  and  UD  aramid  fabrics.  Although  this  study 
showed  that  the  UD  fabric  panels  absorbed  12.5%-  16.5%  more  energy  of  the  projectile  than  did 
the  woven  fabric  for  unit  panel  weights,  the  fibers  used  were  not  the  same. 

In  the  automotive  industry,  Jacob  et  al.16  reviewed  the  energy  absorption  of  polymer 
composites  for  human  safety  in  automotive  crash  events.  It  was  found  that  the  energy  absorption 
by  the  fiber  was  dependent  on  the  fiber  density.  In  addition,  the  energy  absorption  capacity 
increased  with  increasing  interlaminar  fracture  toughness.  As  for  the  fiber-volume  ratio,  as  the 
interlaminar  strength  decreased,  interlaminar  cracks  formed  at  lower  loads,  which  resulted  in  a 
reduced  energy  absorption  capability.  The  density  of  the  composite  increases  with  increasing 
fiber-volume  fraction,  which  results  in  a  lower  energy  absorption  capability  of  the  composite. 

Extensive  research  has  been  conducted  to  address  energy  absorption  behaviors  in  rigid 
UD  composites  subjected  to  various  dynamic  loads;  however,  limited  knowledge  exists 
regarding  the  capacities  of  energy  absorption  mechanisms  in  composites  constructed  with 
compliant  matrix  materials.  Moreover,  the  current  literature  lacks  sufficient  descriptions  of  the 
dominant  failure  mechanisms,  failure  energy  thresholds,  and  limits  for  severe  transient  loadings. 

The  purpose  of  this  research  was  to  demonstrate  that  the  dominant  failure  mechanisms 
that  absorb  the  most  energy  are  dependent  on  the  type  of  dynamic  loads  and,  therefore,  the 
designer  must  fully  know,  a  priori,  the  loading  definitions  for  developing  a  robust  FE  model  for 
predicting  the  dynamic  response  of  the  composite  structure. 


FAILURE  MODES 


The  prevailing  consensus  of  the  scientific  community  considers  fiber  breakage,  matrix 
shearing/cracking,  and  fiber/matrix  debonding  and  delamination  to  be  the  primary  failure 
mechanisms  of  composites.  With  proper  selection  of  the  constituent  materials  and  their 
respective  strengths,  stiffnesses,  and  interfacial  (cohesive)  strengths,  the  designer  can  tailor  the 
thresholds,  limits,  and  sequencing  of  failure  modes  to  enhance  the  dynamic  energy  absorptions 
and  stress  wave  propagations  in  the  composite  structure. 
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Fiber-reinforced  composites  are  often  designed  so  that  the  fibers  bear  the  in-plane  normal 
forces  if  the  composite  is  impacted  in  the  longitudinal  direction.  In  this  scenario,  therefore,  the 
fibers  may  be  first  to  fail.  In  the  event  that  the  fibers  are  loaded  with  some  prescribed  force,  a 
shear  stress  will  develop  down  the  length  of  the  fiber  in  the  interfacial  (cohesive)  zone.  If  the 
shear  stress  induced  from  the  impact  is  greater  than  the  shear  strength  of  the  cohesive  bond,  the 
bond  will  fail  and  will  release  strain  energy — leading  to  a  higher  concentration  of  shear  stress  at 
the  particular  point.  This  fiber/matrix  debonding  is  more  prominent  if  the  bond  is  relatively 
weak.  If,  however,  the  bond  is  relatively  strong,  the  matrix  will  fail  due  to  shear  and  the  failure 
will  extend  along  the  matrix  parallel  to  the  fiber.  As  the  fiber  continues  to  break  in  different 
locations,  the  fiber  will  eventually  delaminate  from  the  matrix  in  its  entirety.  Once  the  fiber 
starts  to  delaminate  from  the  matrix,  the  classic  fiber  pull-out  problem  becomes  an  issue,  as 
described  by  shear  lag  theory.4, 6  Once  the  fiber  is  completely  delaminated  from  the  matrix,  the 
matrix  transfers  the  load  to  the  surrounding  fibers.  There  are  many  theories  describing  the 
process  of  load  transfer  to  neighboring  fibers,  a  review  of  which  can  be  found  in  Brondsted’s4 
article. 


Modeling  the  interfacial  zone  requires  an  understanding  of  cohesive  failure  mechanics 
(that  is,  fiber/matrix  debonding  and  delamination)  and  traction  versus  separation  cohesive 
elements/surfaces.  Cohesive  surfaces  were  used  because  it  was  assumed  that  the  cohesive  layer 
resembles  a  thin  film  with  zero  thickness.  To  implement  cohesive  surfaces  into  the  FE  model,  a 
surface  was  created  around  the  outside  of  the  fiber  and  a  surface  was  created  around  the  inside  of 
the  matrix  where  the  fiber  is  nested.  Initially,  these  two  surfaces  are  bonded  together  and  have  a 
cohesive  stiffness  in  the  normal  ( Kn  )  and  two  shear  directions  ( Ks  and  Kt ).  By  using  a  damage 

initiation  and  evolution  law,  degradation  and  eventual  failure  of  the  bond  between  both  cohesive 
surfaces  can  be  simulated.  Once  the  bond  fails  at  a  particular  interfacial  point,  the  two  surfaces 
at  that  point  are  no  longer  bonded  together  and  cannot  adhere  to  any  other  point  on  either 
cohesive  surface.  In  other  words,  once  an  interfacial  point  fails  (that  is,  the  bond  breaks),  the 
adhesive  characteristics  at  that  point  are  lost  and  irreversible. 

The  contact  stress  or  traction  t  is  defined  as  the  contact  force  divided  by  the  current  area 
at  each  contact  (interface)  point.  The  separation  is  merely  just  the  separation  distance  between 
the  contact  surfaces. 

A  linear  elastic  traction-separation  damage  model  is  used  to  model  the  cohesive  failure 
modes  for  the  cohesive  surfaces  as  shown  in  figure  1 .  The  slope  of  the  first  part  of  the  curve  is 
Kn  (Ks ,  Kt ),  which  represents  the  normal  and  tangential  cohesive  surface  stiffness  components, 
respectively.  Refer  to  table  1  for  the  stiffness  values. 
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Figure  1.  Surface-Based  Cohesive  Behavior  to  Model  the  Interfacial  Zone 


(The  contact  stresses  at  damage  intuition  are  t°n  ,  t°  j  and  effective  separation  at  damage  initiation  8°n  (f> " ,  8°  )■  The 

separations  at  failure  are  Sf  ,  8/  j  where  subscripts  n,  s,  and  t  refer  to  the  normal  and  tangential  directions, 

respectively.) 


The  linear  damage  initiation,  up  to  8°  in  figure  1 ,  is  based  on  a  maximum  separation 
criterion,  which  is  expressed  in  equation  (1)  as 


(1) 


where  8n  is  the  separation  in  the  pure  normal  direction, 

8S  is  the  separation  in  the  first  shear  direction, 

8,  is  the  separation  in  the  second  shear  direction, 

8°,  8%  and  8°  are  the  effective  separations  at  damage  initiation,  and 
l^j  represents  Macaulay  brackets. 

Once  the  separation  criterion  is  met  in  equation  (1)  (refer  to  table  1),  the  contact  point  between 
both  cohesive  surfaces  has  separated  so  that  the  process  of  degradation  begins  by  the  assigned 
damage  evolution  law. 

A  linear  damage  evolution  law  based  on  displacements  is  used  in  this  study  and  is 
represented  by  the  second  half  of  the  curve  in  figure  1 .  As  mentioned,  the  damage  evolution  for 
cohesive  surfaces  describes  the  degradation  of  the  cohesive  stiffness.  A  scalar  damage 
variable/)  represents  the  overall  damage  at  the  contact  point  with  an  initial  value  of  zero.  The 
contact  stress  components  are  a  function  of  the  damage  variable  and  are  defined  in  equations  (2) 
through  (4): 
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otherwise  (no  damage  to  compressive  stiffness) 


If. 

(,=(1-D>“  (2) 

<,=(1  ~Df„ 


where  tn,ts,  and  tt  are  the  contact  stress  components  without  damage.  The  value  of  D  is 


is  the  maximum  value  of  the  effective  separation  attained  during  the  loading  history,  and  S'n  -  8°m 
is  specified  as  the  difference  between  the  separation  at  failure  and  initiation  at  a  contact  point. 


Table  1.  Stiffness  and  Displacement  Properties  for  the  Cohesive  Surfaces 


Kn 

Kn(psi) 

Ks 

Ks(psi) 

Kt 

(psi) 

K 

(in.) 

% 

(in.) 

K 

(in.) 

Sf 

n 

(in.) 

Sf 

n 

(in.) 

8f 

n 

(in.) 

8.0E5 

8.0E5 

8.0E5 

1.7E-5 

1.7E-5 

1.7E-5 

3.0E-5 

3.0E-5 

3.0E-5 

The  same  approach  as  was  described  for  the  cohesive  surfaces  was  taken  to  model  both 
fiber  breakage  and  matrix  shearing — except  for  how  damage  initiation  is  defined.  Damage 
initiation  is  based  on  a  shear  failure  strain  or  an  equivalent  plastic  strain.  The  equivalent  plastic 
strain,  a  scalar  value,  is  defined  in  ABAQUS  as  shown  in  equation  (5): 

t 

s  =  jspdt,  (5) 

o 

where  £p  is  the  plastic  strain  rate  and  defined  in  equation  (6)  as 
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The  plastic  strain  rate  is  deviatoric  with 


£11  =£x2=~>  =  -  |  therefore  2 =  0, 


(V) 


and  all  shear  strain  rates  are  zero.  Once  the  equivalent  plastic  strain  value  is  reached  for  damage 
initiation,  the  same  displacement-based  evolution  law,  as  described  previously  for  the  cohesive 
damage  evolution,  is  used. 


TWO-PLY  QUARTER-SYMMETRY  FINITE  ELEMENT  MODEL 


A  FE  model  of  a  quarter-symmetric  (0/90)  UD  fibrous  composite  laminate  was  developed 
and  simulated  using  ABAQUS/Explicit.17  Two  separate  models,  each  with  different  matrix 
materials,  were  used.  The  first  model  was  a  soft,  compliant  polymer  (Kraton-D),  and  the  second 
was  a  hard  epoxy  resin  (refer  to  tables  2  and  3  for  material  properties).  Each  model  was 
reinforced  with  UHMWPE  fibers  whose  material  properties  are  listed  in  table  4. 


Table  2.  Hyperelastic  Kraton-D  Matrix  Material  Properties 


E 

P  , 

Ultimate  Strength 

ep 

Sf 

E(psi) 

(lb/in.3) 

(psi) 

(in./in.) 

(in.) 

See  figure  2 

8.59E-5 

2800 

1.5 

2.0E-5 

Table  3.  Isotropic  Material  Properties  for  UHMWPE  Fibers18 


E 

o 

P 

Ultimate  Strength 

8f 

(psi) 

(lb/in.3) 

(psi) 

(in./in.) 

(in.) 

17,110,000.0 

.3 

1.3E-4 

580,000 

0.36 

1.0E-5 

Table  4.  Epoxy  Resin  Matrix  Material  Properties 


E 

o 

P 

Ultimate  Strength 

Sp 

8f 

(psi) 

(lb/in.) 

(psi) 

(in./in.) 

(in) 

435,113.0 

A 

1.68E-4 

13,900 

2.5 

2.0E-5 
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The  Kraton-D  polymer  is  a  styrene-isoprene-styrene  (SIS)  hyperelastic  elastomer.  The 
stress  elongation  curve  is  depicted  in  figure  2;  its  ultimate  strength  of  2800  psi  at  1200% 
elongation  was  included  in  the  material  model  for  matrix  failure.  The  hyperelastic  strain  energy 
potential  was  obtained  using  ABAQUS  by  a  third-order  Marlow  model.  The  epoxy  resin  was 
modeled  as  an  isotropic  material  with  an  ultimate  strength  of  13,900  psi.  The  UHMWPE  fibers, 
25  per  ply,  were  given  a  diameter  of  0.001 1  inch  (28  pm)  and  a  length  of  0.025  inch  (635  pm); 
they  were  modeled  as  an  isotropic  material.  The  authors  acknowledge  that  these  fibers  exhibit 
orthotropic  behavior;  however,  the  failure  mechanisms  in  ABAQUS  were  limited  to  isotropic 
material  behavior. 

All  elements  for  both  the  fiber  and  matrix  were  modeled  with  eight-node  linear  brick 
continuum  elements  with  incompatible  modes  (ABAQUS  type  C3D8I).  The  characteristic 
length  of  an  element  is  0.0001  inch. 

The  models  shown  in  the  present  analyses  were  constructed  using  quarter-symmetric 
representations  that  were  permissible  provided  that  the  material  properties,  displacement 
boundary  conditions,  and  applied  loads  were  collectively  quarter- symmetric.  In  figure  3,  the 
quarter- symmetry  planes  are  the  YZ-  and  YX-planes.  The  boundary  conditions  were  such  that 
no  translations  normal  to  the  planes  were  admissible.  The  key  advantage  to  using 
quarter- symmetric  representations  is  reduced  computational  time  with  no  loss  of  solution 
robustness.  Because  of  the  relatively  small  dimensions  of  the  fiber  cross  sections,  ply 
thicknesses,  and  the  level  of  mesh  discretization  applied,  the  solution  time  increment  was  on  the 
order  of  4E-1 1  second. 

The  ABAQUS/Explicit  solver  computes  a  solution  time  increment  based  on  the 
dilatational  wave  speed  in  the  materials  used  and  the  smallest  characteristic  element  length 
according  to 


A  f « 


L 


c 


d 


where  Lmm  is  the  characteristic  length  of  the  smallest  element  and  cd  is  the  dilatational  wave 
speed  defined  as 


where p  is  the  material  density  and  E  is  Young’s  modulus. 

General  contact  surfaces  and  surface  interactions  (hard  overpressure)  were  established 
between  the  matrix  and  fibers  and  between  the  fibers  with  each  other.  For  the  ballistic  impact 
case,  contact  was  defined  for  the  matrix  and  fragment  simulating  projectile  (FSP)  as  well  as  the 
fibers  and  FSP.  It  was  necessary  to  establish  contact  between  the  fibers  and  FSP  because,  as  the 
matrix  material  degraded  due  to  failure,  the  FSP  came  into  contact  with  the  fibers.  If  fiber-to- 
FSP  contact  was  not  established,  the  FSP  would  simply  pass  through  the  fibers. 
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Figure  2.  Hyperelastic  Stress  Elongation  Curve  for  the  Kraton-D  Matrix 


Figure  3.  Left:  Surface  View  of  the  Two-Ply  Model  of  the  Kraton-D  Matrix  (green)  and  the 
UHMWPE  Fibers  (red);  Right:  Subset  View  of  the  Model  Depicting 
the  Mesh  Detail  Adjacent  to  Two  Fibers 
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DYNAMIC  LOADING:  SHOCK 


The  shock  event  was  described  by  applying  an  initial  velocity  to  the  model  at  the  nodes 
along  the  YZ-plane  in  the  negative  X-direction.  The  magnitude  of  the  initial  velocity  was  24,000 
in./s  (2000  ft/s). 


DYNAMIC  LOADING:  BLAST 

The  blast  event  was  described  by  applying  an  overpressure  to  the  upper  lamina  along  the 
XZ-plane  in  the  negative  Y-direction.  The  overpressure  time -history  is  depicted  in  the  upper 
plot  of  figure  4  and  defined  by  equation  (8).  The  equivalent  force-time  history  (based  on 
XZ-surface  area)  is  shown  in  the  lower  plot  of  figure  4.  The  maximum  peak  pressure  was 
3000  psi,  which  decayed  exponentially  for  0.5  ps  where  it  reaches  a  sustained  value  of  107.0  psi 
for  the  remaining  duration  of  the  solution. 

P(f)  =  3000exp(-f/1.5).  (8) 

Because  the  model  was  relatively  small,  achieving  any  deformation  required  a  significant 
pressure  to  be  used  considering  that  the  surface  area  was  0.000625  in.2 — resulting  in  an  applied 
force  of  only  1.875  lb. 


Time  (ns) 


Figure  4.  Upper  Plot:  Overpressure  Curve  for  Dynamic  Blast  Loading  Where  the  Peak 
Pressure  Was  3000  psi  and  Decayed  Exponentially  for  0.5  ps  to  107.0  psi;  Lower  Plot: 
Equivalent  Force  for  the  Blast  Overpressure  Event 


10 


DYNAMIC  LOADING:  BALLISTIC  IMPACT 


To  simulate  a  ballistic  impact  event,  a  quarter-symmetric  representation  of  a  non- 
deformable  FSP  was  used.  The  FSP  was  composed  of  steel  and  was  modeled  as  a  rigid  body 
using  six-node  linear  triangular  prism  continuum  elements  (ABAQUS  type  C3D6).  The  material 
properties  for  steel  are  listed  in  table  5.  A  reference  node  on  the  FSP  was  driven  by  an  initial 
velocity  of  12,000  in./s  (1000  ft/s).  The  assumed  static  coefficients  of  friction  used  between  the 
fibers,  matrix,  and  FSP  was  ps  =0.1. 


Table  5.  Material  Properties  for  the  Steel  FSP 


E  (psi) 

p  (lb/in.3) 

o 

Mass  (lb) 

3.0E7 

7.686E-4 

1.0E-3 

1.313E-10 

RESULTS 


The  FE  models  described  in  this  report  are  representative  of  soft  and  hard  UD  composite 
laminates.  The  matrix  materials  consisted  of  a  hard  epoxy  and  soft  Kraton-D,  and  the  fibers 
were  restricted  to  UHMWPE  fibers.  Both  laminate  models  were  subjected  to  identical  blast, 
shock,  and  ballistic  loading  events  to  investigate  their  energy  absorption  capabilities  and  damage 
mechanisms.  Four  failure  mechanisms  including  fiber  breakage,  matrix  shearing,  and 
fiber/matrix  interfacial  debonding  and  delamination  were  considered.  Each  model  was  simulated 
for  1.2  ps. 

Figures  5  through  7  depict  the  Von  Mises  stress  contours  for  the  blast,  shock,  and 
ballistic  loading  cases  for  the  epoxy  (hard)  matrix,  respectively;  figures  8  through  10  depict  the 
Von  Mises  stress  contours  for  the  blast,  shock,  and  ballistic  loading  cases  for  the  Kraton-D  (soft) 
matrix,  respectively.  For  these  six  illustrations,  the  plots  on  the  right  show  only  the  fiber 
elements,  and  the  plots  on  the  left  show  only  the  matrix  elements. 

The  Von  Mises  stress  contours  for  the  epoxy  matrix  models  blast  and  ballistic  impact 
cases  (see  figures  5  and  7),  achieved  maximum  values  nearly  equal  to  the  specified  ultimate 
strengths  for  both  the  matrix  and  fiber  materials.  For  the  fibers,  the  ultimate  strength  was 
reached  in  the  ballistic  case  (figure  7)  and  is  nearly  developed  in  the  blast  case  (figure  6).  The 
peak  stress  level  in  the  shock  case  did  not  achieve  the  ultimate  strength  of  the  matrix  and  fiber. 

The  Von  Mises  stress  in  the  Kraton-D  matrix  did  not  reach  the  ultimate  strength  in  the 
loading  cases  that  were  considered  because  of  the  material’s  hyperelasticity;  specifically,  the 
hyperelasticity  of  the  Kraton-D  matrix  required  that  the  material  strain  to  1200%  to  achieve  its 
ultimate  stress.  The  ultimate  stress  of  the  fibers  was  reached  in  both  the  blast  (figure  8)  and 
ballistic  (figure  10)  impact  cases.  The  maximum  Von  Mises  stress  of  the  fibers  during  the  shock 
event  (figure  9)  nearly  achieved  the  ultimate  stress. 


11 


S,  Mises 
(Avg:  75%) 

+1 ,390e+04 
+1.275e+04 
+1.159e+04 
+1 ,043e+04 
+9.279e+03 
+8.124e+03 
+6.969e+03 
+5.813e+03 
+4.658e+03 
+3.503e+03 
+2.347e+03 
+1.192e+03 
+3.681e+01 


Figure  5.  Von  Mises  Stress  Contours  in  the  Epoxy  Matrix  (left)  and  UHMWPE  Fibers  (right) 
for  the  Blast  Loading  Case  (Shown  at  Time  t  =  1.2  ps) 


S.  Mises 
(Avg:  75%) 

S+1.128e+04 
+1 ,034e+04 
+9.406e+03 
+8.467e+03 
+7.528e+03 
+6.590e+03 
+5.651  e+03 
+4.713e+03 
+3.774e+03 
+2.835e+03 
+1.897e+03 
+9.583e+02 
+1.974e+01 


S.  Mises 
(Avg:  75%) 

B+2.914e+04 
+2.673e+04 
+2.432e+04 
+2.191e+04 
+1 .950e+04 
+1 ,708e+04 
+1.467e+04 
+1.226e+04 
+9.846e+03 
+7.433e+03 
+5.021  e+03 
+2.609e+03 
+1.961e+02 


Figure  6.  Von  Mises  Stress  Contours  in  the  Epoxy  Matrix  (left)  and  UHMWPE  Fibers  (right) 
for  the  Shock  Loading  Case  (Shown  at  Time  t- 1.2  ps) 
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S,  Mises 
(Avg:  75%) 

E+1 .390e+04 
+1.275e+04 
+1.160e+04 
+1 ,045e+04 
+9.304e+03 
+8.155e+03 
+7.005e+03 
+5.856e+03 
44.707e+03 
+3.558e+03 
+2.409e+03 
+1.260e+03 
+1.106e+02 


S,  Mises 
(Avg:  75%) 

E+5.802e+05 
+5.319e+05 
+4.836e+05 
+4.353e+05 
+3.870e+05 
+3.387e+05 
+2.904e+05 
+2.421e+05 
+1.938e+05 
+1 .455e+05 
+9.726e+04 
+4.897e+04 
+6.789e+02 


Figure  7.  Von  Mises  Stress  Contours  in  the  Epoxy  Matrix  (left)  and  UHMWPE  Fibers  (right) 
for  the  Ballistic  Loading  Case  (Shown  at  Time  t  =  1.2  jus) 


S,  Mises 
(Avg:  75%) 

E+3.333e+03 
+3.055e+03 
+2.777e+03 
+2.500e+03 
+2.222e+03 
+1.944e+03 
+1.666e+03 
+1.389e+03 
+1.111e+03 
+8.333e+02 
+5.555e+02 
+2.778e+02 
+8.962e-02 


S,  Mises 
(Avg:  75%) 

S+5.801e+05 
+5.318e+05 
44.835e+05 
44.353e+05 
+3.870e+05 
+3.387e+05 
+2.904e+05 
+2.421e+05 
+1.939e+05 
+1 ,456e+05 
+9.729e+04 
44.900e+04 
+7.222e+02 


Figure  8.  Von  Mises  Stress  Contours  in  the  Kraton-D  Matrix  (left)  and  UHMWPE  Fibers 
(right)  for  the  Blast  Loading  Case  (Shown  at  Time  t  =  1.2  jus) 
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S,  Mises 
(Avg:  75%) 

■  +1.914e+03 
+1 ,754e+03 
+1.595e+03 
+1.435e+03 
—1-  +1.276e+03 

- +1.116e+03 

—  -  +9.568e+02 

- +7.974e+02 

—  -  +6.379e+02 
-  +4.784e+02 

■  +3.190e+02 
+1.595e+02 
+3.548e-02 


S.  Mises 
(Avg:  75%) 

(+5.140e+04 
+4.712e+04 
+4.284e+04 
+3.856e+04 
+3.428e+04 
+2.999e+04 
+2.571  e+04 
+2.143e+04 
+1.715e+04 
+1.287e+04 
+8.592e+03 
+4.31 1e+03 
+3.060e+01 


Figure  9.  Von  Mises  Stress  Contours  in  the  Kraton-D  Matrix  (left)  and  UHMWPE  Fibers 
(right)  for  the  Shock  Loading  Case  (Shown  at  Time  t  =  1.2  jus) 


S.  Mises 
(Avg:  75%) 

S+9  348e+03 
+8  569e+03 
+7  790e+03 
+7.011e+03 
+6  232e+03 
+5  453e+03 
+4  674e+03 
+3  895e+03 
+3.1 16e+03 
+2.337e+03 
+1  558e+03 
+7.791e+02 
+1  540e-01 


S,  Mises 
(Avg:  75%) 

B+5  802e+05 
+5  318e+05 
+4.835e+05 
+4.352e+05 
+3.868e+05 
+3.385e+05 
+2  901e+05 
+2  418e+05 
+1  935e+05 
+1  451e+05 
+9.681  e+04 
+4.847e+04 
+1  379e+02 


Figure  10.  Von  Mises  Stress  Contours  in  the  Kraton-D  Matrix  (left)  and  UHMWPE  Fibers 
(right)  for  the  Ballistic  Loading  Case  (Shown  at  Time  t  =  1.2  jus) 
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Figure  1 1  depicts  the  energy  dissipated  by  damage  for  the  epoxy  matrix  (A  -  C)  and  the 
Kraton  matrix  (D  -  F)  for  the  blast  (top  (A  and  D)),  shock  (middle  (B  and  E)),  and  ballistic 
(bottom  (C  and  F))  loading  events.  The  inset  plot  in  (A)  shows  the  slight  increase  in  the  matrix 
shearing  (blue  line)  curve.  The  four  failure  mechanisms  are  depicted  as  fiber  breakage  (green 
line),  matrix  shearing  (blue  line),  and  fiber/matrix  cohesive  failure  (red  line).  Cohesive  failure 
was  a  dominant  energy  dissipation  mechanism  for  both  the  blast  and  shock  cases.  There  was  no 
energy  dissipation  from  fiber  breakage  or  matrix  shearing  in  the  shock  case.  Minimal  energy 
dissipation  through  matrix  shearing  occurred  in  the  blast  case,  and  it  was  only  evident  near  the 
end  of  the  simulation  (blue  line). 


Figure  11.  Energy  Loss  Due  to  Damage  for  the  Epoxy  Matrix  (A  -  C)  and  the 
Kraton-D  Matrix  (D  -  F)  for  the  Blast  (top),  Shock  (middle), 
and  Ballistic  Impact  (bottom)  Loading  Events 
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All  failure  modes  contributed  to  the  energy  dissipation  for  the  ballistic  case,  with  fiber 
breakage  being  the  most  dominant.  Matrix  shearing  played  a  larger  role  in  dissipating  energy  in 
the  epoxy  matrix  (C)  than  in  the  Kraton-D  matrix  (F)  for  the  ballistic  case  because  of  the 
elastomeric  behavior  of  the  Kraton-D  matrix,  which  must  severely  elongate  before  failure 
develops.  The  cohesive  failure  dissipates  roughly  the  same  amounts  of  energy.  To  illustrate  this 
further,  refer  to  figures  5  through  9,  which  illustrate  the  composite  laminate  following  the 
dynamic  loading  event.  Note  that,  in  both  the  blast  (figures  5  and  8)  and  shock  (figures  6  and  9), 
very  little  fiber  or  matrix  failure  was  observed — a  finding  also  made  apparent  when  the 
developed  Von  Mises  stress  was  less  than  the  ultimate  strength  of  the  matrix.  For  the  ballistic 
case,  in  both  the  epoxy  matrix  (figure  7)  and  Kraton-D  matrix  (figure  10),  the  projectile  plunged 
through  the  composite  laminate.  All  the  fiber  and  matrix  elements  failed  in  the  area  of  the 
projectile.  In  addition,  the  fibers  delaminated  from  the  matrix  along  the  boundary  of  the 
composite.  These  two  factors  are  why  the  energy  dissipation  due  to  damage  in  the  ballistic 
impact  case  was  much  higher  than  those  produced  in  the  blast  and  shock  cases. 

Figure  12  depicts  the  cohesive  surfaces  for  the  top  ply  (top  row)  and  bottom  ply  (bottom 
row)  for  the  fiber/matrix  interfacial  debonding  and  delamination  for  the  epoxy  matrix  material 
for  blast  (column  A),  shock  (column  B),  and  ballistic  impact  (column  C).  Likewise,  figure  13 
depicts  the  cohesive  surface  damage  results  for  the  Kraton-D  matrix  material.  In  these  figures,  a 
value  of  0  (blue)  indicates  that  the  threshold  values  for  debonding  have  not  been  reached,  and, 
therefore,  cohesive  failure  had  not  yet  intitiated.  A  value  of  1  (red)  indicates  that  a  fiber/matrix 
contact  pair  has  completely  delaminated. 

In  both  matrix  materials  for  the  blast  (column  A)  and  ballistics  (column  C)  case,  the  fiber 
and  matrix  have  delaminated  from  each  other  at  almost  all  interfacial  points.  There  are  slighly 
more  cohesive  failures  in  the  blast  and  ballistic  impact  cases  for  the  Kraton-D  matrix  (figure  13) 
than  for  the  epoxy  case  (figure  12)  because  of  the  high  elongation  to  failure  afforded  by  the 
Kraton-D  matrix.  The  cohesive  failures  for  the  shock  case  (column  B,  figures  12  and  13)  were 
more  uniform  across  the  length  of  the  composite  with  the  Kraton-D  matrix,  but  more  failures 
were  observed  in  the  epoxy  case  for  the  top  ply.  Additionally,  for  the  shock  case,  there  were 
more  cohesive  failures  in  the  bottom  ply  of  the  Kraton-D  matrix  (column  B,  bottom,  figure  13) 
than  the  bottom  ply  of  the  epoxy  matrix  (bottom,  figure  12). 
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Figure  12.  Depiction  of  the  Cohesive  Surfaces  for  the  Top  Ply  (top  row)  and  Bottom  Ply 
(bottom  row)  for  the  Epoxy  Material  Subjected  to  Blast  (column  A), 

Shock  (column  B),  and  Ballistic  Impact  (column  C) 

(The  scale  ranges  from  0-1  where  a  value  of  1  means  that  the  contact  pair  between  that  fiber  and 
matrix  have  delaminated  from  each  other.  A  value  ofO  means  that  the  contact  pair  has  not  yet 
experienced  the  onset  of  damage). 
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Figure  13.  Depiction  of  the  Cohesive  Surfaces  for  the  Top  Ply  (top  row)  and  Bottom  Ply 
(bottom  row)  for  the  Kraton-D  Matrix  Material  Subjected  to  Blast  (column  A), 
Shock  (column  B),  and  Ballistic  Impact  (column  C) 

(The  scale  ranges  from  0-1  where  a  value  of  1  means  that  the  contact  pair  between  that  fiber  and 
matrix  have  delaminated  from  each  other.  A  value  ofO  means  that  the  contact  pair  has  not  yet 
experienced  the  onset  of  damage). 
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Figure  14  depicts  the  velocity-time  history  of  the  projectile  for  the  Kraton-D  (blue  line) 
and  epoxy  matrices  (red  line),  respectively.  The  projectile  in  both  models  penetrated  through  the 
composite  laminate  at  approximately  0.7  ps.  The  projectile’s  residual  velocity  was  less  for  the 
epoxy  matrix  model  versus  the  Kraton-D  matrix  model — suggesting  that  the  epoxy  matrix- 
reinforced  laminate  achieved  a  greater  level  of  energy  absorption. 


Figure  14.  Projectile  Velocity  of  the  Kraton-D  (blue  line)  and  Epoxy  (red  line)  Matrix 


Figure  15  depicts  the  percentage  of  the  total  energy  absorbed  by  the  composite  laminate 
from  the  projectile  impact  for  the  epoxy  (red  line)  and  Kraton-D  (blue  line)  matrices.  The  total 
energy  absorbed  by  the  laminates  consisted  of  the  sum  of  the  elastic  strain  energy,  plastic  strain 
energy,  kinetic  energy,  damage  dissipation  energy,  frictional  dissipation  energy,  and  viscous 
dissipation  energy.  In  agreement  with  figure  14,  the  epoxy  matrix  composite  absorbed  slightly 
more  energy  than  did  the  composite  with  the  Kraton-D  matrix.  The  two-ply  epoxy  matrix  model 
absorbed  4.22%  more  of  the  energy  in  the  same  amount  of  time  as  did  the  two-ply  model  with 
the  Kraton-D  matrix  for  this  particular  ballistic  impact  event. 

Time -history  plots  of  each  active  energy  absorption  mechanism  are  presented  in 
figure  16.  Figure  16  (top)  depicts  the  individual  energy  absorption  terms  for  both  the  epoxy 
(column  A)  and  Kraton-D  (column  B)  matrix  subjected  to  ballistic  loading.  The  energy  terms 
include  damage  dissipation  (blue  line),  frictional  dissipation  (green  line),  kinetic  energy  of  the 
composite  (red  line),  elastic  strain  energy  (magenta  line),  plastic  strain  energy  (cyan  line),  and 
viscous  dissipation  (yellow  line).  The  frictional  energy  was  based  on  the  use  of  static  friction 
coefficients  defined  for  contact  between  the  projectile,  fibers,  and  matrices.  The  viscous 
dissipation  utilized  a  default  value  that  ABAQUS  included. 
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Figure  15.  Percentage  of  the  Projectile  Kinetic  Energy  Absorbed  by  the  Epoxy  (red  line) 
and  Kraton-D  (blue  line)  Composite  Laminate 
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Figure  16.  Top  Row:  Percentage  of  the  Energy  Absorption  Mechanisms  of  the  Composite 
Laminate  Including  Damage  Dissipation  (blue  line),  Frictional  Dissipation  (green  line), 
Kinetic  Energy  of  the  Composite  (red  line),  Plastic  Dissipation  (cyan  line),  Strain  Energy 
(magenta  line),  and  Viscous  Dissipation  (yellow  line)  for  the  Epoxy  (column  A)  and  Kraton-D 
(column  B)  Matrix  Materials  Subjected  to  Ballistic  Loading;  Bottom  Row:  Energy  Balance 
(cyan  line)  Showing  the  Sum  of  the  Energy  Absorbed  by  the  Composite  (blue  line)  and  the 

Kinetic  Energy  of  Projectile  (red  line) 
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The  sum  of  all  energy  terms  is  represented  by  the  blue  line  in  figure  16  (bottom);  the 
kinetic  energy  of  the  projectile  (red  line),  a  replicate  of  figure  15,  was  added  here  for 
convenience.  The  sum  of  these  two  plots  yields  a  constant  value  (cyan  line)  over  the  course  of 
the  simulation.  Because  this  energy  curve  is  constant,  the  time  rate  of  change  would  yield  a 
value  of  zero,  indicating  conservation  of  energy.  Effectively,  the  figure  is  graphically  showing 
that  equation  (9)  is  being  satisfied,  where  the  blue  line  is  the  right  side  of  the  equation  and  the 
difference  of  the  cyan  line  and  red  line  is  the  left  side  of  the  equation. 
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In  both  matrix  composite  configurations,  the  dictating  energy  absorption  terms  are  the  kinetic 
energy,  elastic  strain  energy,  and  plastic  strain  energy  of  the  composite. 


DISCUSSION 


The  purpose  of  this  research  was  to  investigate  the  energy  absorption  capabilities  of 
admissible  failure  modes  for  a  cross-ply  UD  composite  laminate  subjected  to  severe  dynamic 
loading  events  including  blast,  shock,  and  ballistic  impact.  The  composite  laminates  used  in  this 
research  were  composed  of  hard  (epoxy)  and  soft  (Kraton-D)  matrix  materials.  The  fibers  used 
in  both  matrix  models  were  UHMWPE  fibers.  The  simulation  time  of  each  loading  case  was 
1.2  ps. 


Damage  in  the  fibers,  matrix,  and  cohesive  interfaces  was  observed  to  be  interdependent 
as  the  load  path  from  the  fiber  to  fiber-matrix  interface  to  matrix  was  sequential.  If  the  fiber  and 
or  matrix  have  failed  at  a  particular  region,  then  cohesive  damage  cannot  initiate  and  accumulate 
at  that  region.  The  same  is  not  true,  however,  for  the  reverse  case.  If  a  cohesive  interface  has 
delaminated,  the  fiber  and  the  matrix  at  that  region  can  still  support  stress  and  possibly  lead  to 
failure;  therefore,  all  the  cohesive  interfaces  can  fail  while  the  fiber  and  matrix  continue  to 
function  and  later  experience  damage  initiation  and  evolution.  This  possibility  is  illustrated  in 
the  epoxy  matrix  case  for  the  blast  loading  event  (refer  to  figure  1 1  (A)  and  figure  12  (column 
A)).  At  the  end  of  the  simulation,  almost  all  the  cohesive  surfaces  have  failed,  which  is  indicated 
by  a  value  of  1  (red).  In  addition,  the  dissipation  due  to  cohesive  damage  plateaued  in  figure  1 1 
(A  (red  line));  however,  at  approximately  1  ps,  with  nearly  all  cohesive  regions  failed,  matrix 
failure  just  initiated — depicted  by  the  small  curve  representing  matrix  shearing  in  figure  1 1  (A 
(blue  line)).  Further  blast  loading  increased  matrix  failures  despite  a  majority  of  the  cohesive 
contact  regions  delaminating. 

Compared  to  fiber  breakage  and  matrix  shearing,  cohesive  damage  was  the  dominant 
damage  mechanism  for  energy  dissipation  in  both  the  blast  and  shock  cases  for  both  the  epoxy 
and  Kraton-D  matrix  material  models.  The  amounts  of  energy  dissipation  through  cohesive 
damage  in  the  epoxy  and  Kraton-D  matrix  models  were  roughly  the  same  for  the  respective  blast 
cases;  however,  the  Kraton-D  matrix  models  achieved  this  amount  at  a  much  faster  rate  than  did 
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the  epoxy  matrix  models.  The  faster  rate  was  attributed  to  the  lower  cohesive  strength  of  the 
Kraton-D  matrix  in  contrast  to  that  of  epoxy.  Fiber- to-matrix  delamination  develops  faster  in 
time  for  matrices  having  lower  cohesive  strengths  in  a  given  blast  event.  This  important 
observation  is  critical  to  the  design  of  protective  composites  for  blast  loading  in  which  fiber-to- 
matrix  cohesive  damage  is  a  necessary  mechanism  for  energy  dissipation. 

Results  of  the  shock  loading  case  revealed  that  the  accumulation  of  cohesive  damage 
occurred  earlier  in  time  and  to  a  greater  extent  for  the  epoxy  matrix  model  than  for  the  Kraton-D 
matrix  model.  Here,  the  fibers  were  loaded  along  their  longitudinal  axis  compressing  the  matrix 
material.  A  matrix  material  that  is  weaker  than  the  fiber  tries  to  peel  itself  back  down  the  length 
of  the  fiber.  Because  the  epoxy  matrix  material  is  more  compressible  than  is  the  elastomeric 
Kraton-D  material,  it  will  experience  more  cohesive  failure  than  will  Kraton-D  material. 

The  dominant  damage  mechanism  for  the  ballistic  impact  case  was  fiber  breakage 
(followed  by  appreciable  matrix  shearing  for  the  epoxy  matrix  model).  As  the  projectile  sheared 
through  the  material  (referred  to  as  “plug  formations”),  the  primary  fibers  and  matrix  elements  at 
the  projectile  contact  region  failed  by  shearing  through  the  laminate.  Furthermore,  additional 
failures  of  the  primary  fibers  on  the  strike  face  and  several  secondary  fibers  on  the  rear  face  were 
observed  along  the  model  boundaries,  thus  increasing  the  total  damage  developed  through  fiber 
breakage.  As  in  the  blast  and  shock  cases,  the  energy  dissipation  due  to  matrix  shearing  damage 
was  notable  in  the  epoxy  matrix  model  and  nearly  negligible  for  the  Kraton-D  matrix  model. 

The  Kraton-D  matrix  did  not  appreciably  resist  shearing  deformations  because  of  its  relatively 
high  compliance  and  low  shear  strength.  With  the  epoxy  matrix,  however,  energy  dissipation 
through  matrix  shearing  damage  was  relatively  low.  Likewise,  the  cohesive  damage  energies  for 
the  Kraton-D  and  epoxy  materials  were  also  nearly  negligible. 

Figures  12  and  13  show  that,  for  the  ballistic  impact  case,  both  matrix  materials 
experienced  cohesive  damage  initiation  nearly  throughout  the  volumes  of  their  respective 
models.  Furthermore,  nearly  all  the  cohesive  interface  regions  experienced  damage  initiation, 
which,  for  the  Kraton-D  matrix,  led  to  full  cohesive  failures.  The  amount  of  cohesive  damage 
energies,  however,  was  nearly  negligible — not  surprising  in  that  there  were  far  fewer  cohesive 
interface  regions  to  fail  compared  to  the  number  of  fiber  and  matrix  elements  that  could  fail  for  a 
given  fiber  volume  ratio  Vf.  Each  fiber  element,  matrix  element,  and  cohesive  region  that  fails 
dissipates  energy  through  damage.  It  is  surmised,  therefore,  that  laminated  composites 
constructed  of  increasing  values  of  Vf  have  a  greater  propensity  than  do  laminated  composites 
constructed  of  decreasing  values  of  V/to  develop  greater  cohesive  damage  energies. 

The  epoxy  matrix  (red  line  in  figure  14)  absorbed  more  ballistic  impact  energy  than  did 
the  Kraton-D  matrix  (blue  line  in  figure  14)  as  evidenced  by  the  projectile’s  smaller  residual 
velocity  Vr,  which  is  2.9%  less.  Further  evidence  that  the  epoxy  matrix  absorbed  more  ballistic 
impact  energy  than  did  the  Kraton-D  matrix  is  shown  in  figure  15,  which  depicts  the  total  kinetic 
energy  of  the  projectile  absorbed  by  the  composite  laminate.  The  energy  absorbed  using  the 
epoxy  matrix  (red  line)  is  4.22%  greater  than  that  absorbed  due  to  damage  by  the  Kraton-D 
matrix  (blue  line). 
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The  upper  time  history  plots  of  figure  16  depict  each  of  the  energy  terms  of  the  ballistic 
impact  case  for  the  epoxy  matrix  (column  A)  and  the  Kraton-D  matrix  (column  B).  The 
dominating  energy  absorption  mechanisms  of  the  composite  were  the  plastic  strain  energy, 
elastic  strain  energy,  and  kinetic  energy.  The  remaining  energy  terms,  roughly  one  order  of 
magnitude  less,  were  the  viscous  energy,  damage  energy,  and  frictional  energy.  Because  energy 
was  conserved,  the  total  sum  of  the  absorbed  energies  (blue  line)  and  the  kinetic  energy  of  the 
projectile  (red  line)  were  plotted  by  a  constant  value  curve  (cyan  line)  in  figure  16  (bottom). 

This  constant  value  was  the  projectile’s  initial  kinetic  energy,  thus  confirming  energy 
conservation. 

Lastly,  the  epoxy  matrix  absorbed  49.17%  of  the  initial  projectile  kinetic  energy  (4.22% 
more  energy  than  that  absorbed  by  the  Kraton-D  matrix);  however,  more  matrix  material  was 
observed  to  have  failed  for  the  epoxy  matrix  than  for  the  Kraton-D  matrix.  This  finding  suggests 
that  full  projectile  arrest  could  be  achieved  by  doubling  the  number  of  plies,  resulting  in  an 
approximate  %-ballistic  limit  and  the  choice  of  matrix  material  will  also  dictate  the  amount  of 
energy  absorbed  which  becomes  pivotal  in  composite  design  and  optimization.  Additionally,  by 
accounting  for  all  energy  absorption  and  failure  mechanisms,  a  manufacturer  could  better 
produce  composite  laminates  by  tailoring  the  matrix  and  fiber  materials;  tensile,  shearing,  and 
cohesive  failure  thresholds;  and  ply  counts  and  stacking  schedules  to  optimally  defeat  specific 
projectile  threats. 


CONCLUSIONS 


This  research  centered  on  investigating  the  energy  absorption  capabilities  of  failure 
modes  in  unidirectional  laminated  composites  subjected  to  various  dynamic  loading  conditions 
including  blast,  shock,  and  ballistic  impact.  Four  failure  mechanisms  were  considered:  fiber 
breakage,  matrix  shearing,  fiber  delamination,  and  matrix  delamination  (that  is,  cohesive  failure). 
Two  matrix  materials,  a  rigid  epoxy  and  a  compliant  elastomer  (Kraton-D),  were  evaluated.  The 
reinforcing  fiber  material  was  ultrahigh  molecular  weight  polyethylene. 

The  results  suggest  that  the  energy  dissipated  through  cohesive  failure  mechanisms  must 
be  accounted  for  in  the  blast  and  shock  events.  The  rate  at  which  energy  was  dissipated  through 
cohesive  damage  for  the  blast  case  differed  for  the  two  matrix  materials;  however,  the  amount  of 
energy  dissipated  through  cohesive  failures  was  nearly  invariant  for  the  two  matrices.  For  the 
shock  event,  the  epoxy  matrix  dissipated  more  energy  through  cohesive  damage  than  did  the 
Kraton-D  matrix.  All  damage  mechanisms  must  be  accounted  for  in  ballistic  impact  events; 
however,  fiber  breakage  was  the  most  dominant  damage  energy  mechanism.  It  was  shown  that 
the  epoxy  matrix  better  minimized  the  projectile’s  residual  velocity  and  absorbed  4.22%  more 
dynamic  energy.  More  matrix  material,  however,  has  failed. 

While  the  threshold  values  of  damage  initiation  and  evolution  used  throughout  this  study 
were  limited  to  single  values,  the  inclusion  of  cohesive  failure  modes  was  shown  to  be  critical  in 
developing  robust  FE  models  for  both  blast  and  shock  events.  Material  properties  including 
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failure  thresholds  must  be  obtained  through  experiments  conducted  on  the  fibers,  matrix,  and  the 
composite  laminate  system. 


The  authors  acknowledge  that  this  research  considered  only  one  fiber  volume  ratio  V/and 
that  future  studies  should  address  the  effects  of  matrix-rich  and  matrix-starved  constructions. 
Similarly,  the  ballistic  impact  case  assumed  a  rigid  (nondeforming)  projectile  with  a  right 
circular  cylinder  shape.  Other  projectile  types  and  velocities  should  be  included  in  future 
studies. 
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